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Abstract 

The hard-core model has attracted much attention across several disciplines, 
representing lattice gases in statistical physics and independent sets in discrete 
mathematics and computer science. On finite graphs, we are given a parameter 
A, and an independent set I arises with probability proportional to A^L On 
infinite graphs a Gibbs distribution is defined as a suitable limit with the correct 
conditional probabilities. In the infinite setting we are interested in determining 
when this limit is unique and when there is phase coexistence, i.e., existence 
of multiple Gibbs states. On finite graphs we are interested in determining the 
mixing time of local Markov chains. 

On Z 2 it is conjectured that these problems are related and that both undergo 
a phase transition at some critical point A c th 3.79 [2]. For the question of phase 
coexistence much of the work to date has focused on the regime of uniqueness, 
with the best result to date being recent work of Restrepo et al. [24] showing 
that there is a unique Gibbs state for all A < 2.3882. Here we give the first 
non-trivial result in the other direction, showing that there are multiple Gibbs 
states for all A > 5.3646. Our proof begins along the lines of the standard Peierls 
argument, but we add two significant innovations. First, building on the idea of 
fault lines introduced by Randall [23], we construct an event that distinguishes 
two boundary conditions and yet always has long contours associated with it, 
obviating the need to accurately enumerate short contours. Second, we obtain 
vastly improved bounds on the number of contours by relating them to a new 
class of self-avoiding walks on an oriented version of Z 2 . 

The best result for rapid mixing of local Markov chains on boxes of Z 2 is 
also when A < 2.3882 [24]. Here we extend our characterization of fault lines 
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to show that local Markov chains will mix slowly when A > 5.3646 on lattice 
regions with periodic (toroidal) boundary conditions and when A > 7.1031 with 
non-periodic (free) boundary conditions. The arguments here rely on a careful 
analysis that relates contours to taxi walks and represent a sevenfold improvement 
to the previously best known values of A [23] . 

1 Introduction 

The hard-core model was introduced in statistical physics as a model for lattice gases, 
where each molecule occupies non-trivial space in the lattice, requiring occupied sites 
to be non-adjacent. When a lattice such as 7L d is viewed as a graph, the allowed 
configurations of molecules naturally correspond to independent sets in the graph. 

Given a graph G, let Q be the set of independent sets of G. Given a (fixed) activity 
(or fugacity) A G M + , the weight associated with each independent set / is w(I) = A' 7 '. 
The associated Gibbs (or Boltzmann) distribution \i = hg,\ is denned on Q, assuming G 
is finite, as //(/) = w(I)/Z, where the normalizing constant Z = Z(G, A) = J2j e nw(J) 
is commonly called the partition function. Physicists are interested in the behavior of 
models on an infinite graph (such as the integer lattice Z d ), where the Gibbs measure 
is denned as a certain weak limit with appropriate conditional probabilities. For many 
models it is believed that as a parameter of the system is varied - such as the inverse 
temperature j3 for the Ising model or the activity A for the hard-core model - the system 
undergoes a phase transition at a critical point. 

For the classical Ising model, Onsager, in seminal work [20], established the precise 
value of the critical temperature (3 C (Z 2 ) to be log(l + \/2). Only recently have the 
analogous values for the (more general) g-state Potts model been established in break- 
through work by Beffara and Duminil-Copin [3], settling a more than half-a-century old 
open problem. Establishing such a precise value for the hard-core model with currently 
available methods seems nearly impossible. Even the existence of such a (unique) criti- 
cal activity A c , where there is a transition from a unique Gibbs state to the coexistence 
of multiple Gibbs states remains conjectural for Z d (d > 2; it is folklore that there is 
no such transition for d — 1), while it is simply untrue for general graphs (even general 
trees, in fact, thanks to a result of Brightwell et al. [7]). Regardless, a non-rigorous 
prediction from the statistical physics literature [2] suggests A c ~ 3.796 for Z 2 . 

Thus, from a statistical physics or probability point of view, understanding the 
precise dependence on A for the existence of unique or multiple hard-core Gibbs states 
is a natural and challenging problem. Moreover, breakthrough works of Weitz [29] 
and Sly [27] in recent years identified A c (Ta) - the critical activity for the hard-core 
model on an infinite A-regular tree - as a computational threshold where estimating the 
hard-core partition function on general A-regular graphs undergoes a transition from 
being in P to being iVP-hard (specifically, there is no PTAS unless NP = RP), further 
motivating the study of such (theoretical) physical transitions and their computational 
implications. While it is not surprising that for many fundamental problems computing 
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the partition function exactly is intractable, it is remarkable that even approximating 
the partition function of the hard-core model above a certain critical threshold also 
turns out to be hard. 

Starting with Dobrushin [8] in 1968, physicists have been developing techniques to 
characterize the regimes on either side of A c for the hard-core model. Most attention has 
focused on establishing ever larger values of A below which there is always uniqueness 
of phase. The problem has proved to be a fruitful one for the blending of ideas from 
physics, discrete probability and theoretical computer science, with improvements to 
our understanding of the problem having been made successively by Radulescu and 
Styer [22], van den Berg and Steif [4] and Weitz [29], among others. The state of the 
art is recent work of Restrepo et al. [24] , building on the novel arguments introduced 
by Weitz, which establishes uniqueness for all A < 2.3882. 

Much less is known about the regime of phase coexistence. Dobrushin [8] established 
phase coexistence for all A > C, but did not explicitly calculate C. Around the time 
of the writing of [6], Borgs reported that C = 80 was the theoretical lower limit of 
Dobrushin's argument [5], but a recent computation by the second author suggested 
that the actual consequence of the argument was more like C ~ 300. 

Our first main result: Using insights from physics and combinatorics, we establish 
the first non-trivial upper bound on A c for the hard-core model on Z 2 . 

Theorem 1.1 For all A > 5.3646, the hard-core model on Z 2 with activity A admits 
multiple Gibbs states. 

From a computational standpoint, there are two natural questions to ask concerning 
the hard-core model on a finite graph. Can the partition function be approximated? 
(an exact determination is in general too much to hope for), and how easy is it to 
sample from a given Gibbs distribution? For both questions, a natural and powerful 
method is offered by Markov chain algorithms - carefully constructed random walks on 
the space of independent sets of a graph whose equilibrium distributions are the desired 
Gibbs distributions. One of the most commonly studied families of Markov chains are 
the local-update chains, such as Glauber dynamics, that change the state of a bounded 
number of vertices at each step (in particular, Glauber dynamics changes the state of 
a single vertex at each step). 

The efficiency of the Markov chain method relies on the underlying chain being 
rapidly mixing; that is, it must fairly quickly reach a distribution that is quantifiably 
close to stationary. For many problems, local chains such as Glauber dynamics seem to 
mix rapidly below some critical point, while mixing slowly above that point. Most no- 
tably for the Ising model on Z 2 , simple local Markov chains are rapidly mixing (in fact, 
with optimal rate) for j3 < (3 C (Z 2 ) and slowly mixing for (3 > /3 C (Z 2 ); various mathe- 
matical physics experts (including Aizenmann, Holley, Martinelli, Olivieri, Schonmann, 
Stroock and Zegarlinski) have contributed to this work. Recently the Ising picture was 
completed by Lubezky and Sly [16], who showed polynomial mixing at j3 — /3 C (Z 2 ). 
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Once again, the known bounds are less sharp for the hard-core model. Luby and 
Vigoda [17] showed that Glauber dynamics on independent sets is fast when A < 1 
on the 2-dimensional lattice and torus. Weitz [29] reduced the analysis on the grid to 
the tree, thus establishing that in this same setting Glauber dynamics is fast up to the 
critical point for the 4-regular tree, in effect for A < 1.6875. The state of the art in this 
direction is due to Restrepo et al. [24]. Using now standard machinery (establishing 
so-called strong spatial mixing), they proved that the natural Glauber dynamics on the 
space of hard-core configurations on boxes in Z 2 is rapidly mixing for all A < 2.3882. 
These results also lead to efficient deterministic algorithms for approximating the hard- 
core partition function on (finite regions of) Z 2 . 

As with phase-coexistence, it is believed that there is a critical value A" ux across 
which the Glauber dynamics for sampling from hard-core configurations on a n by n 
box in Z 2 flips from mixing in time polynomial in n, to exponential in n, and that 
it coincides with A c . Pinning down such a precise transition seems beyond reach at 
the moment, but we can try to establish a short range of values into which A™ lx (if it 
exists) must fall, by finding ever smaller bounds for A above which the mixing time is 
exponential. 

Borgs et al. [6] showed that Glauber dynamics is slow on toroidal lattice regions in Z d 
(for d > 2), when A is sufficiently large, in particular establishing a finite constant above 
which mixing is slow on Z 2 . The first effective bound was provided by Randall [23], 
who showed slow mixing for A > 50.59 on boxes with periodical boundary conditions, 
and for A > 56.88 on boxes with free boundary. (Better bounds were originally reported 
in [23], but these were due to a minor error; in any case, the bounds we report in the 
present paper improve on the original claims as well.) 

Our second main result: Here we apply the combinatorial insights gained from the 
proof of Theorem 1.1 to establish slow mixing of Glauber dynamics on boxes in Z 2 
for values of A that are an order of magnitude lower than the previously best known 
bounds. 

Theorem 1.2 For all A > 5.3646, the mixing time for Glauber dynamics for the hard- 
core model on n by n boxes in Z 2 with periodic boundary conditions and with activity 
A is exponential in n. For free boundary conditions, we have the same result with 
A > 7.1031. 

Our techniques: The standard approach to showing multiple Gibbs distributions in a 
statistical physics model is to consider the limiting distributions corresponding to two 
different boundary conditions on boxes in the lattice centered at the origin, and find a 
statistic that separates these two limits. For the hard-core model, it suffices (see [4]) to 
compare the even boundary condition - all vertices on the boundary of a box at an even 
distance from the origin are occupied - and its counterpart the odd boundary condition, 
and the distinguishing statistic is typically the occupation of the origin. Under odd 
boundary condition the origin should be unlikely to be occupied, since independent 
sets with odd boundary and (even) origin occupied must have a contour - a two-layer 
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thick unoccupied loop of vertices separating an inner region around the origin that is in 
"even phase" from an outer region near the boundary that is in "odd phase" . For large 
enough A, such an unoccupied layer is costly, and so such configurations are unlikely. 
This is essentially the Peierls argument for phase coexistence, and was the approach 
taken by Dobrushin [8]. 

As we will see presently, the effectiveness of the Peierls argument is driven by the 
number of contours of each possible length - better upper bounds on the number of 
contours translate directly to better upper bounds on A c . Previous (unpublished) work 
on phase coexistence in the hard-core model on Z 2 had viewed contours as simple 
polygons in Z 2 , which are closely related to the very well studied family of self-avoiding 
walks. While this is essentially the best possible point of view when applying the Peierls 
argument on the Ising model, it is far from optimal for the hard-core model. One of the 
two major breakthroughs of the present paper is the discovery that hard-core contours, 
if appropriately defined, can be viewed as simple polygons in the oriented Manhattan 
lattice (orient edges of Z 2 that are parallel to the x-axis (resp. y-axis) positively if 
their y-coordinate (resp. ^-coordinate) is even, and negatively otherwise), with the 
additional constraint that contours cannot make two consecutive turns. The number of 
such polygons can be understood by analyzing a new class of self-avoiding walks, that 
we refer to as taxi walks. The number of taxi walks turns out to be significantly smaller 
than the number of ordinary self-avoiding walks, leading to much better bounds on A c 
than could possibly have been obtained previously. 

There is a single number fi t > 0, the taxi walk connective constant, that asymptot- 
ically controls the number c n of taxi walks of length n, in the sense that c n = ji^f t {n) 
with ft(n) sub-exponential. Adapting methods of Aim [1] we obtain good estimates 
on fi t , giving us good understanding of c n for large n. The sub-exponential correction 
makes it difficult to control c n for small n, however, presenting a major stumbling block 
to the effectiveness of the Peierls argument. Using the statistic "occupation of origin" 
to distinguish the two boundary conditions, one inevitably has to control c n for both 
small and large n. The lack of precise information about the number of short contours 
leads to discrepancies between theoretical lower limits and actual bounds, such as that 
between C = 80 (theoretical best possible) and C « 300 for phase coexistence on Z 2 
discussed earlier. 

The second breakthrough of the present paper is the idea of using an event to 
distinguish the two boundary conditions that has the property that every independent 
set in the event has associated with it a long contour. This allows us to focus exclusively 
on the asymptotic growth rate of contours/taxi walks, and obviates the need for an 
analysis of short contours. The immediate result of this breakthrough is that the 
actual limits of our arguments agree exactly with their theoretical counterparts. The 
distinguishing event we use extends the idea of fault lines, which we discuss in more 
detail below in the context of slow mixing. 

The traditional argument for slow mixing is based on the observation that when A 
is large, the Gibbs distribution favors dense configurations, and Glauber dynamics will 
take exponential time to converge to equilibrium. The slow convergence arises because 
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the Gibbs distribution is bimodal: dense configurations lie predominantly on either the 
odd or the even sublattice, while configurations that are roughly half odd and half even 
have much smaller probability. Since Glauber dynamics changes the relative numbers 
of even and odd vertices by at most 1 in each step, the Markov chain has a bottleneck 
leading to torpid (slow) mixing. 

Our present work builds on a novel idea from [23] in which the notion of fault 
lines was introduced to establish slow mixing for the Glauber dynamics on hard-core 
configurations for moderately large A, still improving upon what was best known at 
that time. Randall [23] gave an improvement by realizing that the state space could 
be partitioned according to certain topological obstructions in the configurations, rather 
than the relative numbers of odd or even vertices. Not only does this approach give 
better bounds on A, but it also greatly simplifies the calculations. First consider an 
n x n lattice region G with free (non-periodic) boundary conditions. A configuration / 
is said to have a fault line if there is a width two path of unoccupied vertices in / from 
the top of G to the bottom or from the left boundary of G to the right. Configurations 
that do not have a fault line must have a cross of occupied vertices in either the even 
or the odd sublattices forming a connected path in G 2 from both the top to the bottom 
and from the left to the right of G, where G 2 connects vertices at distance 2 in G. 
Roughly speaking the set of configurations that have a fault line forms a cut set that 
must be crossed to move from a configuration that has an odd cross to one with an 
even cross, and it was shown that fault lines are exponentially unlikely when A is large. 
Likewise, if G is an n x n region with periodic boundary conditions, it was shown that 
either there is an odd or an even cross forming non-contractible loops in two different 
directions or there is a pair of non-contractible fault lines, and a similar argument can 
be made. 

We improve this argument by refining our consideration of fault lines. Previous 
bounds just characterized fault lines as (rotated) self-avoiding walks in Z 2 ; here we 
observe that, suitably modified, they are in fact taxi walks and so the machinery de- 
veloped for phase coexistence can be brought to bear in the mixing context, leading to 
our improvements. 

We believe that both of the breakthroughs of the present paper have more general 
applicability. A natural next step is to study the hard-core model on Z d for d > 3. 
Establishing reasonable bounds on the critical activity for Z 3 is an immediate challenge, 
as is pinning down how the criticaljvalue changes with d. The best upper bounds are 
0(<i~ 1//4 ) for slow mixing [10] and 0(rf -1 / 3 ) [21] for phase coexistence; the best known 
lower bounds in both cases are Q(d~ r ). 

The rest of this manuscript is laid out as follows. Section 2 provides much of the 
necessary combinatorial background material. This includes the notion of fault lines 
and crosses, the connection between fault lines and taxi walks, and the characterization 
of independent sets on finite regions of Z 2 . In Section 3 we explain how fault lines and 
taxi walks can be used to prove Theorem 1.1 (phase coexistence). In Section 4 we prove 
Theorem 1.2 (slow mixing). Finally, in Section 5 we give the detailed analysis of the 
growth rate of taxi walks. 
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2 Combinatorial background 



Here we introduce the notions of crosses, fault lines and taxi walks, which are the key 
ingredients of the proofs of both Theorem 1.1 and Theorem 1.2. 

2.1 Crosses and fault lines 

We begin by defining some useful graph structures. Let G = (V, E) be a simply con- 
nected region in Z 2 , say the nxn square. We define the graph Go = (V<>, E$) as follows. 
The vertices V$ are associated with the midpoints of edges in E. Vertices u and v in Vq 
are connected by an edge in E§ if and only if they are the midpoints of incident edges 
in E that are perpendicular. Notice that Go is a region in a smaller Cartesian lattice 
that has been rotated by 45 degrees. 

We will also make use of the even and odd subgraphs of G. For b G {0, 1}, let 
Gf, = (Vb, Eb) be the graph whose vertex set Vj, C V contains all vertices with parity 
b (i.e., the sum of their coordinates has parity b), with (u,v) G Ef, if u and v are 
connected in G 2 . We refer to Go and G\ as the even and odd subgraphs. The graphs 
Go, Go and Gi play a central role in defining the features of independent sets that 
determine distinguishing events in our proof of phase coexistence and the partition of 
the state space for our proofs of slow mixing. 

Given an independent set I G Q, we say that a simple path p in Go is spanning if it 
extends from the top boundary of Go to the bottom, or from the left boundary to the 
right, and each vertex in p corresponds to an edge in G such that both endpoints are 
unoccupied in I. It will be convenient to color the vertices in V$ along a spanning path 
using the parity of the vertex to the "left" (or "top") of the path in V. In particular, 
recall that each vertex v G V$ on the path p bisects an edge e v G E. Each edge in 
E has an odd and an even endpoint, and we color v blue if the odd vertex in e v is 
to the left when the path crosses v, and red otherwise. Every time the color of the 
vertices along the path changes, we have an alternation point. It was shown in [23] 
that if an independent set has a spanning path, then it must also have one with zero 
or one alternation points. We call this path a fault line, and we let fij- be the set of 
independent sets in Q that contain at least one fault line. 

We say that I G has an even bridge if there is a path from the left to the right 
boundary or from the top to the bottom boundary in Go consisting of occupied vertices 
in 1. Similarly, we say it has an odd bridge if it traverses Gi in either direction. We say 
that I has a cross if it has both left-right and a top-bottom bridges. 

Notice that if an independent set has an even top-bottom bridge it cannot have an 
odd left-right bridge, so if it has a cross, both of its bridges must have the same parity. 
We let f2 Q ^ be the set of configurations that contain an even cross and let fii C f2 
be the set of those with an odd cross. 

We can now partition the state space Q into three sets, with one separating the 
other two; this partition is critical to the proofs of both Theorem 1.1 and Theorem 1.2. 
The following lemma was proven in [23]. 
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Lemma 2.1 Let G be as above. The set f2 of independent sets on G can be partitioned 
into sets Q^, Qo and Qi, consisting of configurations with a fault line, an even cross or 
an odd cross. If I <E Q and I' G Q± then \IAI'\ > 1. 

It will be useful to extend these definitions to the torus as well. Let n be even, and 
let G be the nx n toroidal region {0, . . . , n — 1} x {0, . . . n — 1}, where v = (vi,v 2 ) and 
u = (ui, u 2 ) are connected if V\ = U\ ± l(mod n) and v 2 = u 2 or v 2 = u 2 ±l(mod n) and 
v 1 — u\. Let Q be the set of independent sets on G and let n be the Gibbs distribution. 
As before, we consider Glauber dynamics that connect configurations with symmetric 
difference of size one. 

We define G<>, Go and G\ as above to represent the graph connecting the midpoints of 
perpendicular edges (including the boundary edges), and the odd and even subgraphs. 
As with G, all of these have toroidal boundary conditions. 

Given J e O, we say that / has a fault F = (F 1: F 2 ) if there are a pair of vertex- 
disjoint non-contractible cycles F 1 ,F 2 in G§ whose vertices correspond to edges in G 
where both endpoints are unoccupied, and such that the vertices on each cycle are all 
red or all blue (i.e., the endpoints in G to one side of either cycle all have the same 
parity). We say that / has a cross if it has at least two non-contractible cycles of 
occupied sites in I with different winding numbers. The next lemma (from [23]) utilizes 
faults to partition Q. 

Lemma 2.2 Let G be as above. The set VI of independent sets on G can be partitioned 
into sets Qjr,Q ,Q 1; consisting of configurations with a fault, an even cross or an odd 
cross. If I G Q and V G Qi then \IAI'\ > I. 

2.2 Taxi walks 

The strategy for the proofs of phase coexistence and slow mixing will be to use a Peierls 
argument to define a map from Qjr to Q that takes configurations with fault lines to 
ones with exponentially larger weight. The map is not injective, however, so we need 
to be careful about how large the pre-image of a configuration can be, and for this it is 
necessary to get a good bound on the number of fault lines. In [23] the number of fault 
lines was bounded by the number of self- avoiding walks in Gq (or G$ on the torus). 
However, this is a gross overcount because, as we shall see, this includes all spanning 
paths with an arbitrary number of alternation points. As we shall see, we can get much 
better bounds on the number of fault lines by only counting a subset of self-avoiding 
walks with zero or one alternation points. 

To begin formalizing this idea, we put an orientation on the edges of Gq. Each 
edge (u, v) G E§ corresponds to two edges in E that share a vertex w G V. We orient 
the edge "clockwise" around w if w is even and "counterclockwise" around w if w is 
odd. For paths with zero alternation points, all of the edges must be oriented in the 
same direction (with respect to this edge orientation). If we rotate so that the 
edges are axis aligned, then this simply means that the horizontal (resp. vertical) edges 
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alternate direction according to the parity of the y- (resp. x-) coordinates, like in many 
well-known metropolises. 

We can now define taxi walks. Let Z 2 be an orientation of Z 2 in which an edge par- 
allel to the x-axis (resp. y-axis) is oriented in the positive x-direction if its ^/-coordinate 
is even (resp. oriented in the positive y-direction if its ^-coordinate is even), and is 
oriented in the negative direction otherwise (note that this corresponds exactly to the 
orientation placed on G§ above). It is common to refer to Z 2 as the Manhattan lattice: 
streets are horizontal, with even numbered streets oriented East and odd numbered 
streets oriented West, and avenues are vertical, with even numbered avenues oriented 
North and odd numbered avenues oriented South. 

Definition 2.1 A taxi walk is an oriented walk in Z 2 that begins at the origin, never 
revisits a vertex, and never takes two left or two right turns in a row. 

We call these taxi walks because the violation of either restriction during a real taxi 
ride would cause suspicion among savvy passengers. 

Lemma 2.3 If an independent set I has a fault line F with no alternations, then it 
also has a fault line F' so that either F' or F' R (the reversal of F') is a taxi walk. 

Proof: It is straightforward to see that if I has a fault line F with no alternation 
points, then it must have all of its edges oriented the same way (in Gq) and it must be 
self-avoiding. Suppose F is a minimal length fault line in I without any alternations, 
and suppose that F has two successive turns. Because of the parity constraints, the 
vertices immediately before and after these two turns must both connect edges that are 
in the same direction, and these five edges can be replaced by a single edge to form a 
shorter fault line without any alternations. This is a contradiction to F being minimal, 
completing the proof. □ 

The same argument can be used to show that if F is a fault line with an alternation 
point, then there is a fault line that is the concatenation of two taxi walks (or the rever- 
sals of taxi walks). Lemma 2.3, and the extension just mentioned, are key ingredients 
in our proofs of both phase coexistence and slow mixing. They allow us to assume, as 
we do throughout, that all fault lines we work with are essentially taxi walks. For phase 
coexistence we will also need to understand the connection between Peierls contours 
and taxi walks. 

Given an independent set I in Z 2 , let (I°) + be the set of odd vertices in I together 
with their neighbors. Let R be any finite component of (I°) + , and let W be the unique 
infinite component of Z 2 \ R. Let C be the complement of W (the process of going from 
R to C is essentially one of "filling in holes" in R). Finally, let 7 be the set of edges 
with one end in W and one in C, and write 70 for the subgraph of induced by 7. 

Lemma 2.4 In G§, 7^ is a directed cycle that does not take two consecutive turns. 
Consequently, if an edge is removed from 7, the resulting path in Z 2 (suitably translated 
and rotated so that it starts at the origin) is a taxi walk. 
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Proof: Because 7 separates W from its complement, 7^ must include a cycle surround- 
ing a vertex of W, and since 7 is in fact a minimal edge cutset (W and C are both 
connected), 7^ must consist of just this cycle. To see both that 7^ is correctly (i.e. 
uniformly) oriented in G<>, and that it does not take two consecutive turns, note that 
if either of these conditions were violated then we must have one of the following: a 
vertex of W (or C) all of whose neighbors are in C (or W), or a unit square in Z 2 with 
both even vertices in C and both odd vertices in W (this is an easy case analysis). All 
of these situations lead to a 4-cycle in 7^, a contradiction since 7^ is a cycle whose 
length is evidently greater than 4 (in fact it must have length at least 12). □ 

A critical step in our arguments will be bounding the number of taxi walks. We 
start by recalling facts about standard self-avoiding walks which have been studied 
extensively, although many basic questions remain (see, e.g., [18]). It is easy to see that 
on Z 2 , the number c n of walks of length n grows exponentially with n as 2™ < c n < 
4 x 3 n_1 , since there are always at most 3 ways to extend a self-avoiding walk of length 
n — 1 and walks that only take steps to the right or up can always be extended in 2 
ways. Hammersley and Welsh [14] showed that c n = \i n exp(0(y/n)), where u is known 
as the connective constant. It is believed (and there is considerable experimental and 
heuristical evidence to suggest) that exp(0(\/n)) here can be replaced by 0(n n / 32 ). 

Let c n be the number of taxi walks of length n. It is easy to see that 2™/ 2 <c n < 
4 x 3 n_1 . For the upper bound here we use c n < c„, and for the lower bound we observe 
that if we take two steps at a time in one direction we can always go East or North. 
With very little extra work, we can get a significantly better upper bound: 

Lemma 2.5 Letc n be the number of taxi walks of length n. Thenc n = 0((1 + v / 5)/2)™. 

Proof: At each vertex there are exactly two outgoing edges in Z 2 . If we arrive at v 
from u, then one of the outgoing edges continues the walk in the same direction and 
the other is a turn. The two allowable directions are determined by the parity of the 
coordinates of v, so we can encode each walk as a bitstring s e {0, l} n_1 . If s± — 
then the walk starts by going East (along a street) and if So — 1 the walk starts North 
along an avenue. For all % > 1, if Sj = the walk continues in the same direction as 
the previous step, while if Sj = 1 then the walk turns in the permissible direction. In 
this encoding, the condition forbidding consecutive turns forces s to avoid having two 
l's in a row, and hence c n < f n = O(0 n ), where f n is the nth Fibonacci number and 
= (1 + V5)/2 m 1.618 is the golden ratio. □ 

Using more sophisticated tools, described in more detail in Section 5, we can improve 
our bounds. 

Theorem 2.6 Let c n be the number of taxi walks of length n. Then c n = 0(1.5884)™. 
Moreover, there exists a constant /i t > ( the taxi walk connective constant ) such that 
for all n, c n = f t (n)n™, where ft{n) grows sub- exponentially. We have 1.5196 < fi t < 
1.5884 and 4.3332 < fij — 1 < 5.3646. If /i > ii t then for all large n, c n < u n . 
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3 Proof of Theorem 1.1 



Here we establish phase coexistence for the hard-core model on Z 2 for all sufficiently 
large A. The following stronger statement implies Theorem 1.1 via Theorem 2.6. 

Theorem 3.1 The hard-core model on Z 2 with activity A admits multiple Gibbs states 
for all A > fif — 1 , where fi t is the connective constant of taxi walks. 

We will not review the theory of Gibbs states, contenting ourselves with saying infor- 
mally that an interpretation of the existence of multiple Gibbs states is that the local 
behavior of a randomly chosen independent set in a box can be made to depend on a 
boundary condition imposed on the box, even in the limit as the size of the box grows 
to infinity. See e.g. [12] for a very general treatment, or [4] for a treatment specific to 
the hard-core model on the lattice. 

Let U n be the box [— n, +n} 2 , and I Q the independent set consisting of all even 
vertices of Z 2 . Let J* be the set of independent sets that agree with J e off U n , and 
jji e n the distribution supported on in which each set is selected with probability 
proportional to A' /nf/n l Define /x° analogously (with "even" everywhere replaced by 
"odd"). We will exhibit an event A that depends only on finitely many vertices, with 
the property that for all large n, n%[A) < 1/3 and [i° a (A) > 2/3. This is well known 
(see e.g. [4]) to be enough to establish the existence of multiple Gibbs states. 

The event A depends on a parameter m = m(X) whose value will be specified later. 
Specifically, A consists of all independent sets in Z 2 whose restriction to U m contains 
either an odd cross or a fault line. We will show that ^{A) < 1/3 for all sufficiently 
large n; reversing the roles of odd and even throughout, the same argument gives that 
under /x° the probability of U m having either an even cross or a fault line is also at most 
1/3, so that (by Lemma 2.1) fi° dd (A) > 2/3. 

Write A e n for A D note that for all large n we have ^(A) = ^(A^). To show 
H e n (A c n ) < 1/3 we will use the fact that / G A c n is in even phase (predominantly even- 
occupied) outside U n , but because of either the odd cross or the fault line in U m it is 
not in even phase close to U m ; so there must be a contour marking the furthest extent 
of the even phase inside U n . We will modify I inside the contour via a weight-increasing 
map, showing that an odd cross or fault line is unlikely. 

3.1 The contour and its properties 

Fix I G A c n . If / has an odd cross in U m , we proceed as follows (using the notation from 
the discussion preceding Lemma 2.4). Let R be the component of that includes 

a particular odd cross. Note that because I agrees with I e off U n , R does not reach the 
boundary of U n , and so as in the discussion preceding Lemma 2.4, we can associate to 
R a cutset 7 separating it from the boundary of U n . 

Notice that 7 is an edge cutset in U n separating an interior connected region that 
meets U m from an exterior connected region that includes the boundary of U n , with all 
edges from the interior of 7 to the exterior going from an unoccupied even vertex to an 
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unoccupied odd vertex. This implies that |7|, the number of edges in 7, is a multiple 
of 4, specifically four times the difference between the number of even and odd vertices 
in the interior of 7. Because the interior includes two points of the odd cross that are 
at distance at least 2m + 1 from each other in U m , we have a lower bound on 7 that is 
linear in m; in particular, clearly I7I > m. Note also that by Lemma 2.4, 7^ is a closed 
taxi walk. 

We now come to the heart of the Peierls argument. If we modify / by shifting it 
by one axis-parallel unit (positively or negatively) in the interior of 7 and leaving it 
unchanged elsewhere, then the resulting set is still independent, and we may augment it 
with any vertex in the interior whose neighbor in the direction opposite to the shift is in 
the exterior. This is a straightforward verification; see [6, Lemma 6] or [11, Proposition 
2.12] where this is proved in essentially the same setting. Furthermore, from [6, Lemma 
5] each of the four possible shift directions free up exactly |7|/4 vertices that can be 
added to the modified independent set. 

We now describe the contour if / has a fault line in U m . If there happens to be an 
odd occupied vertex in U m then we construct 7 as before, starting with some arbitrary 
component of (/ ) + that meets U m in place of the component of an odd cross. If the 
resulting 7 has a fault line in its interior, then 7 and its associated 7^ satisfy all the 
previously established properties immediately. 

Otherwise, choose a fault line, which we can assume by Lemma 2.3 is a taxi walk 
or the concatenation of two taxi walks. Whether it has zero or one alternation points, 
we can find a path P = UiUi . . . Uk in Z 2 with k linear in m, with U\ and both odd, 
with no two consecutive edges parallel, and with the midpoints of the edges of the path 
inducing an alternation-free sub-path of the chosen fault line (essentially we are just 
taking a long piece of the fault line, on an appropriately chosen side of the alternation 
point, if there is one). This sub-path F 1 is a taxi walk. Next, we find a second path 
in Gq, disjoint from F ± , that always bisects completely unoccupied edges, and that 
taken together with Fi completely encloses P. If there are no occupied odd vertices 
adjacent to even vertices of P, such a path is easy to find: we can shift F\ one unit in 
an appropriate direction, and close off with an additional edge at each end. If there are 
some odd occupied vertices adjacent to some even vertices of P, then this translate of 
Fi has to be looped around the corresponding components of Such a looping is 

possible because does not reach the boundary of U n , nor does it enclose the fault 

line (if it did, we would be in the case of the previous paragraph). 

This second path we have constructed may not be a taxi walk; however, following 
the proof of Lemma 2.3, we see that a minimal path F 2 satisfying the conditions of our 
constructed path is indeed a taxi walk. We take the concatenation of i*\ and F 2 to be 
7^ in this case, and take 7 to be the set of edges that are bisected by vertices of 7^. 
The contours in this case satisfy all the properties of those in the previous case. The 
standard strategies outlined in [6] and [1 1] can easily be used to derive the properties in 
this case. The one difference is that now 7^ may not be a closed taxi walk; but at worst 
it is the concatenation of two taxi walks, both of length linear in m (and certainly it 
can be arranged that each has length at least m/2). 
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3.2 The Peierls argument 

For J G J* set w(J) = A |JnC H Our aim is to show that w(A e n )/w(J°) < 1/3. For 
/ G A e n) let ip(I) be the set of independent sets obtained from / by shifting in the 
interior parallel to (1,0) and adding all subsets of the |7|/4 vertices by which the 
shifted independent set can be augmented. For J G let S denote the set of added 
vertices. Define a bipartite graph on partite sets A e n and J* by joining / G A c n to 
J G J* if J G </?(/). Give edge JJ weight u>(7)A 1 ' 5 ' 1 = iu(J) (where 5 is the set of 
vertices added to / to obtain J). 

The sum of the weights of edges out of those I G A c n with |7(/)| = 4£ is (1 + A)^ 
times the sum of the weights of those I. For each J G J*, the sum of the weights of 
edges into J from this set of Ps is w(J) times the degree of J to the set. If f(£) is a 
uniform upper bound on this degree, then 



The lower bound on i here is crucial. The standard Peierls argument takes A to be the 
event that a fixed vertex is occupied, and the analysis of probabilities associated with 
this event requires dealing with short contours, leading to much weaker bounds than 
we are able to obtain. 

To control f{£), observe that for each J G and contour 7 of length 4£ there is 
at most one / with 7(f) =7 such that J G (p(I) (I can be reconstructed from J and 
7, since the set S of added vertices can easily be identified; cf. [11, Section 2.5]). It 
follows that we may bound f(£) by the number of contours of length M with a vertex 
of U m in their interiors. 

Fix fi > fif By the properties of contours we have established, up to translations of 
contours this number is at most the maximum of fx M and Ylj+k=At j k>m/2 V^V^ = ^\x u 
(for all large m, here using Theorem 2.6). The restriction of to U m has at most 
4(2m + l) 2 < 17m 2 edges, so there are at most this many translates of any particular 
contour that can have a vertex of U m in its interior. If follows that we may bound 
f(£) by 68m 2 £/x 4 ^ and so the sum in (1) by X^> m 68m 2 £(/x 4 /(l + A))^. For any fixed 
A > fi 4 — 1, there is an m large enough so that this sum is at most 1/3; we take any 
such m to be m(A), completing the proof of phase coexistence. 

4 Proof of Theorem 1.2 

We now bound the mixing time of Glauber dynamics by showing that the conductance 
is exponentially small. 

4.1 Glauber dynamics and mixing time 

Let G C Z 2 be an n x n lattice region and let Q be the set of independent sets on G. 
Our goal is to sample from Q according to the Gibbs distribution, where each I G is 




(1) 
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assigned probability 

tt(J) = \^/Z, 

and Z = X^/'en ' * s the normalizing constant. 

Glauber dynamics is a local Markov chain that connects two independent sets if they 
have symmetric difference of size one. The Metropolis probabilities [19] that force the 
chain to converge to the Gibbs distribution are given by 

f^min(l,Al / 'l-l / l) , if I@I' = 1, 

[0, otherwise. 

The conductance, introduced by Jerrum and Sinclair [26], is a good measure of the 
mixing rate of a chain. Let 

, . 52 xe s,vts*( x ) p ( x >y) 
<p = mm — 

Sen:7r(5)<l/2 7l(S) 

where ir(S) = J2 xeS n(x) * s tne weight of the cutset S. The following classical theorem 
provides the connection between low conductance and slow mixing. 

Theorem 4.1 [26] For any Markov chain with conductance $ we have ^ < Gap(P) < 
2$, where Gap(P) is the spectral gap of the transition matrix. 

The spectral gap is well-known to be a measure of the mixing rate of a Markov chain 
(see, e.g., [25]), so an exponentially small conductance is sufficient to show slow mixing. 
Using Theorem 4.1, our goal is therefore to define a partition that has exponentially 
small conductance. The machinery developed in Section 2.1 (in particular Lemmas 2.1 
and 2.2) provide exactly such a partition. 

4.2 Glauber dynamics on the 2-d torus 

We are now ready to complete the proof of slow mixing, starting first with the two- 
dimensional torus. Let n be even, and let G = {0, . . . , n — 1} x {0, . . . , n — 1} be the 
n x n lattice region with toroidal boundary conditions. We take Q to be the set of 
independent sets on G and let n be the Gibbs distribution. Lemma 2.2 shows that 
fl may be partitioned into three sets, Qjr (independent sets containing a fault), Q 
(independent sets with an even cross) and Q\ (independent sets with an odd cross), 
and that furthermore f2 and Q\ are not directly connected by moves in the chain P. 
It remains to show that 7r(Qjr) is exponentially smaller than both 7r(Sl ) an d 
(Clearly tt(Q q ) = tt(Qi) by symmetry.) Notice that on the torus we may assume that 
fault lines have no alternation points; since they start and end at the same place, the 
number of alternation points must be even. 

Given an independent set / G Qjr with fault F = (F\, F2), we partition the vertices 
of / into two sets, I a and Ib, depending on which side of Fi and F 2 they lie. Define 
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the length of a fault to be the total number of edges on this path (or paths) in G<>. 
Notice that all fault lines with zero alternation points have length N = n + 2£, for some 
integer £, since they all have the same parity. 

Let V = a (I, F) be the configuration formed by shifting I a one to the right. Let 
F[ = cr(Fi) and F 2 = o~(F 2 ) be the images of the fault under this shift. We define the 
points that lie in Fi fl F[ and F 2 fl F 2 to be the points that fall "in between" F and 

It will be convenient to order the set of possible fault lines so that given a config- 
uration I e Qjr we can identify its first fault. The following results are modified from 
[23] and rely on the new characterization of faults as taxi walks. 

Lemma 4.2 Let Vl F be the configurations in fijr with first fault F = (F 1 ,F 2 ). Write 
the length of F as 2n + 2£. Then 

tt(q f ) < (i + \y( n+i \ 

Proof: We define an injection <p F :Vt F x {0, l} n+e ^ so that n((f) F (I,r)) = tt(I)X^. 
The injection is formed by cutting the torus G along Fi and F 2 and shifting one of 
the two connected pieces in any direction by one unit. There will be exactly n + £ 
unoccupied points near F that are guaranteed to have only unoccupied neighbors. We 
add a subset of the vertices in this set to I according to bits that are one in the vector 
r. 

Given this map, we have 

1 = 7?(fi) 

/en F re{o,i}™+ £ 

= £?(/) E ^ H 

i<=h F re{o,i}"+ e 

□ 

Theorem 4.3 Let fi be the set of independent sets on G weighted by = X^/Z, 
where Z = Yliefi Let fijr be the set of independent sets on G with a fault. Then 
for any A > /4 — 1, there is a constant c > such that 

7r(JV) < e" cn . 

Proof: Fix ji satisfying A > /i 4 — 1 > /i 4 — 1, where \i t is the taxi walk connective 
constant. Summing over possible locations for the two faults F\ and F 2 and using 
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Lemma 4.2, 

n(Q T ) = E*W 



< ^(l + A)-(™ +£ ) 

F 

(ri 2 -2n)/2 . . 
,_n V / 



i=0 

4 \ n+i 



2 \^ / A* 



E 



l + A 



The second inequality here uses Theorem 2.6. By our choice of \x we get (for large n) 
^(^V) < e ~ cn f° r some constant c > 0; and we can easily modify this constant to deal 
with all smaller values of n. □ 

From Section 2.2 we know that fif — 1 < 5.3646. Combining Theorems 4.1 and 4.3, 
we thus get the first part of Theorem 1.2, as well as the following stronger result. 

Corollary 4.4 Fix A > fif — 1. Glauber dynamics for sampling independent sets on 
the n x n torus G takes time at least e cn to mix, for some constant c > (depending 
on A ). 

Proof: We will bound the conductance by considering the cut S = VLq. It is clear that 
< 1/2 since S = Qjr U Qi and tt(Qq) = tt(D,i). Thus, 

$ < $ 5 

E fle n ,ten^( s ) P (M) 

7r(^0) 

< E fe sy gg) 

7r(^0) 
ggV) 
TT(^o) ' 

Given Theorem 4.3, it is trivial to show that tt(£1 ) > 1/3, thereby establishing that the 
conductance is exponentially small. It follows from Theorem 4.1 that Glauber dynamics 
takes exponential time to converge. □ 



4.3 Non-periodic boundary conditions 

For regions with non-periodic boundary conditions we also employ a weight-increasing 
map from configurations with fault lines by performing a shift and adding vertices. In 
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this setting, however, we not only have to reconstruct the position of the fault line, but 
we must also encode the part of the configuration lost by the shift due to the finite 
boundary. In this section we give the proof of the following result, which implies the 
second part of Theorem 1.2. 

Theorem 4.5 Fix A satisfying 

1. 1 + A > n\ and 

2. 2(1 + A) > /i 2 t (l + VT+IX), 

where ji t is the taxi walk connective constant. Glauber dynamics for independent sets 
on the n x n grid G takes time at least e cn to mix, for some constant c = c(A) > 0. 

Simple algebra reveals that the second condition is satisfied whenever A 2 + (2 — // 2 — 
/4)A + (1 — > 0- Using fx t < 1.5884 we find that both of these conditions are met 
when A > 7.1031, as claimed in Theorem 1.2. 

As before, we partition the state space Q into three sets, namely Qj- (independent 
sets with a fault line), Qq (independent sets with an even cross), and Qi (independent 
sets with an odd cross). From Lemma 2.1 we know that these sets do indeed partition 
Q, and that furthermore Q$ and are not connected by moves in P. It remains to 
show that vr(fijr) is exponentially smaller than 7r(f2 ) and 7r(Oi) . 

Let I G Qjr be an independent set with vertical fault line F. The fault F partitions 
the vertices of G into two sets, Right (F) and Left(F), depending on the side of the 
fault on which they lie. Recall that a fault has zero or one alternation points, and the 
edges form a path in G<>. We will represent length of the path as N = n + 2£, for some 
£ G Z. Strictly speaking this is only correct for fault lines with zero alternation points; 
but this slight abuse will affect the analysis of what follows by only a constant factor. 

Let V = a (I, F) be the configuration formed by shifting Right (F) one to the right. 
We will not be concerned right now if some vertices "fall off" the right side of the 
region G. Let F' = cr(F) be F shifted one to the right. We define the points that lie in 
Right (F) n Left(F') to be the points that fall "in between" F and F' . 

We use the following result from [23] . 

Lemma 4.6 Let I be an independent set with a fault line F. Let V = a(F,I) and 
F' = cr(F) be defined as above. 

1. F and F' are both fault lines in V . 

2. If we form V 1 by adding all the points that lie in between F and F' to V ( except 
the unique odd point incident to the alternation point, if it exists), then V' will be 
an independent set. 

3. If \F\ = n + 2£, then there are exactly n + £ points that lie in between F and F' . 

Moreover, F and F' are taxi walks or the concatenation of two taxi walks at an alter- 
nation point. 
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Let / G f2jr be an independent set with a fault line, which we assume is vertical. 
(If / only has horizontal fault lines, we can rotate G so that it is vertical; the net 
effect of ignoring these independent sets is at most a factor of 2 in the upper bound on 
7r(f2j-), and this will get incorporated into other constant factors.) Let F = F(I) be 
the leftmost fault line. Let the length of F be n + 2£, for some integer £. 

Let Gi tTl be the lxn lattice representing the last column of G, and let J be any 
independent set on Gi jU . We further partition Qj?, into Uf,j^f,j, where I G &f,j if it 
has leftmost fault line F and is equal to J when restricted to the last column G\ yn . 

Lemma 4.7 Let F be a fault in G with length n + 2£ and let 5 equal the number of 
alternation points on F (so 5 = or 1). Let J be an independent set on Gi tn . With 
^f,j defined as above, we have 

n(n F ,j) < Al J l(l + A)-< n+< -*>. 

Proof: Let r G {0, iy n + e ~ s be any binary vector of length n + £ — 5 and let |r| denote 
the number of bits set to 1, where |r| < n + £. The main step is to define an injective 
map 4> F ,j : ^Ijr x {0, l} n+e — y Q such that, for any I G IV, 

7r(0 F)J (/,r)) = 7r(/)A-l J l + l r l. 

Given this map, we have 

1 = 7r(Q) 

ien FiJ re{o,i} n+e - s 

= E E -(^)a- |j|+m 

ien FtJ re{o,i} n+e - s 

= E ^ lJl E Alr| 

ien FiJ re{o,i} n+e - s 

= vr(/)A- |J| (l + A) n+ ^ 

ien FrJ 

= x~\ J \(i + X) n+i ~ 5 ir(n FtJ ). 

We define the injective map <Pf,j in stages. For any / G &f,j, w e delete the last 
column (which is equal to J). Next, recalling that any fault line partitions G into two 
pieces, we identify all points in / that fall on the right half and shift these to the right by 
one using the map a(I, F). From Lemma 4.6 we know that the number of points that 
fall between these two fault lines is n + £, where n + 2£ is the length of the fault. The 
final step defining the map is to insert new points into the independent set along this 
strip between the two faults using the vector r, thereby adding |r| new points. The new 
independent set (f>F,j(I,r) has |/| — \ J\ + \r\ points, and hence has weight 7r(/)A~' J ' + ' r ' 

□ 
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Lemma 4.8 Let G^ n be a 1 x n strip, and let £l n be the set of independent sets in G^ r . 
Then 



E AW < c (i±4±5y 

IcO. V / 



for some constant c. 

Proof: Let be the set of independent sets in G\ yi and let Tj = ^Jen* ^' J '- Then 
T = 1, 7\ = 1 + A, and 

Tj = Tj_! + ATj_ 2 . 

Solving this Fibonacci-like recurrence yields the lemma. □ 

Theorem 4.9 Let Q be the set of independent sets on the n x n lattice G weighted by 
n (I) = \^/Z, where Z = J2 Ien A' 7 ' is the normalizing constant. Let f2jr be the set of 
independent sets on G with a fault line. Then 

tt(JV) < p(n) e~ c ' n , 

for some polynomial p{n) and constant c' > 0, whenever A satisfies the hypothesis of 
Theorem 4-5. 

Proof: We will make use of the injective map <f)p t j : &f,j x {0, 1}^ — > where 
N = n + 2£ is the length of the fault line. Using Lemma 4.8 for the third inequality 
and Theorem 2.6 for the fourth, we have 

F,J 

< ^Al J l(l + A)-(" + ^ 

F,J 

< A^(l + A)^^Al J l 

F Jefir 

1 + VTT4AA" 



< Ae E (l + A)-"+'>( 1 + V ^y 

2 

i=o V - / 

= AC "?ITTaJ( 2(1 + A) )■ 



This means that we will have 7r(f2jr) < p{n)e c ' n , for some polynomial p(n), if both 

(l + A)> / u 4 , 2(1 + A) > n 2 (l + Vl + 4A) 

hold. Taking arbitrarily close to fit, the result follows. □ 
Theorem 4.5 follows exactly as Corollary 4.4 in Section 4.2. 
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5 Taxi walks: Bounds and Limits 



We conclude by justifying the upper bound on the number of taxi walks given in The- 
orem 2.6, as well as providing a lower bound on ji t . It is necessary to first establish the 
submultiplicativity of c n (or, equivalently, the subadditivity of logc n ). 

Lemma 5.1 Let c n be the number of taxi walks of length n and let 1 < i < n — 1 . Then 

Proof: As with traditional self-avoiding walks, the key is to recognize that if we split 
a taxi walk of length n into two pieces, the resulting pieces are both self-avoiding. Let 
s — si, . . . , s n be a taxi walk of length n and let 1 < % < n — 1. Then the initial segment 
of the walk sj = s\, . . . , Sj+i is a taxi walk of length i. Let p = (x, y) be the ith vertex 
of the walk s. Let sf be the final n — i steps of the walk s starting at p. We define 
f(sp) by translating the walk so that f(p) is the origin, reflecting horizontally if p x is 
odd and reflecting vertically if p y is odd. Notice that this always produces a valid taxi 
walk of length n — i and the map / is invertible given p. Therefore c n < Ci c n ^i. □ 

It follows from Lemma 5.1 that a n = logc n is subadditive, i.e., a n+m < a n + a m . By 
Fekete's Lemma (see, e.g., [28, Lemma 1-2.2]) we know that lim^oo a n /n exists and 

lim^ = inf^. (2) 

n->co n n 

Thus, we can write the number of taxi walks as c n = fi™ft(n), where \i t is the connective 
constant associated with taxi walks and ft{n) is subexponential in n. 

Subadditivity gives us a strategy for getting a better bound on /i t . From (2) we see 
that for all n, \ogc n /n is an upper bound for log// t . We exactly enumerated taxi walks 
of length n, for n < 60; see http://nd.edu/~dgalvinl/TD/ for this and other data. 
Using cqq = 2189670407434 gives a bound of fx t < 1.6058. Note that exact counts for 
larger n will immediately improve our bounds on both fi t and A c . 

The connective constant for ordinary self-avoiding walks has been well studied, and 
some of the methods used to obtain bounds there can be adapted to deal with taxi 
walks. In particular, a method of Aim [1] is useful. Fix n > m > 0. Construct a square 
matrix A(m,n) whose ij entry counts the number of taxi walks of length n that begin 
with the ith taxi walk of length to, and end with the jth taxi walk of length to, for 
some fixed ordering of the walks of length to. To make sense of this, it is necessary to 
choose, for each v e Z , an orientation preserving map /„ of Z that sends the origin 
to v; saying that a walk of length n ends with the jth walk of length to means that if 
the length to terminal segment of the walk is transformed by f^ 1 to start at the origin, 
where v is the first vertex of the terminal segment, then the result is the jth walk of 
length to. Then a theorem of Aim [1] may be modified to show that fi t is bounded 
above by Xi(A(m, n)) 1 ^™ - " 1 ^, where Ai indicates the largest positive eigenvalue. (Note 
that when to = this recovers the subadditivity bound discussed earlier). 

We have calculated A(20,60). This is a square matrix of dimension 20114, and a 
simple symmetry argument reduces the dimension by a factor of 2. Using MATLAB, 
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we could estimate the largest eigenvalue of this reduced matrix to obtain \x t < 1.5884 
and nf-K 5.3646. 

A similar strategy can be used to derive lower bounds on \i t in order to determine 
the theoretical limitations of our approach of characterizing contours by taxi walks. 
We have already given the trivial lower bound (j, t > \/2. To improve this, we consider 
bridges (introduced for ordinary self-avoiding walks by Kesten [15]). A bridge, for our 
purposes, is a taxi walk that begins by moving from the origin (0, 0) to the point (1, 0), 
never revisits the y-axis, and ends by taking a step parallel to the x-axis to a point on 
the walk that has maximum x-coordinate over all points in the walk (but note that this 
maximum does not have to be uniquely achieved at the final point). 

Let b n be the number of bridges of length n. Then bridges are supermultiplicative, 
i.e., b n > bib n -i (and logb n is superadditive). To see this, note that if fa and fa are 
bridges, then they both begin and end at vertices whose y-coordinates are even because 
they are taking steps to the East. If the parities of the x-coordinates of the first 
vertices in fa and fa agree, then the concatenation of fa and an appropriate translation 
of fa is also a bridge; if the parities are different then concatenation of fa with a 
translation of fa after reflecting horizontally will be a valid bridge. Notice that the 
parity of the x-coordinate of the two pieces allows us to recover whether a reflection 
was necessary to keep the walk on the directed Manhattan lattice, so bridges are indeed 
supermultiplicative. It similarly follows that there are at least b k taxi walks of length 
kn (just concatenate k length n bridges), so that 

fit = lim c]l m > lim (b k n ) l/nk = b\[ n . 

m— s-oo A;— »oo 

We have enumerated bridges of length up to 60, in particular discovering that &60 — 
80312795498, leading to fi t > 1.5196 and $ - 1 > 4.3332. 

A consequence of our lower bound on ji t is that our present approach to phase 
coexistence cannot give anything better than A c < 4.3332; this tells us that new ideas 
will be needed to reach the value of 3.796 suggested by computations as the true value 
of A c . 
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